NON UNIQUE SOLUTIONS TO BOUNDARY VALUE PROBLEMS 
FOR NON SYMMETRIC DIVERGENCE FORM EQUATIONS 



ANDREAS AXELSSON 



Abstract. We calculate explicitly solutions to the Dirichlet and Neumann bound- 
ary value problems in the upper half plane, for a family of divergence form equa- 
tions with non symmetric coefficients with a jump discontinuity. It is shown that 
the boundary equation method and the Lax-Milgram method for constructing so- 
lutions may give two different solutions when the coefficients are sufficiently non 
symmetric. 



1. Introduction 

Recently, new techniques in harmonic analysis have been used to study boundary 
value problems (EVP's) for divergence form elliptic equations with non symmetric, 
or more general complex coefficients. In the half plane, for real but non symmet- 
ric coefficients, Lp solvability of the Dirichlet problem for sufficiently large p was 
obtained by Kenig, Koch, Pipher and Toro [3] and Lp solvability of the Neumann 
and regularity problems, for sufficiently small p, was proved by Kenig and Rule [4]. 
In R", two boundary equation methods have been studied by Alfonseca, Auscher, 
Axelsson, Hofmann and Kim f\] and by Auscher, Axelsson and Hofmann [2] where, 
among other things, these EVP's are proved to be well posed in L2 for small complex 
Loo perturbations of real symmetric coefficients. 

Unlike the case of real symmetric coefficients, for general non symmetric coeffi- 
cients, well posedness of these classical EVP's may fail. In [3] and [1], the family 



Ak(x) 



1 ksgn{x) 
—ksgn{x) 1 



of non symmetric coefficient matrices with a jump at x = was studied and shown 
to provide counter examples to well posedness for certain values of the parameter 
/c e R. More precisely, the following theorem was proved in (3], Theorem (3.2.1)] 
and [H Appendix]. 

Theorem 1.1. Let 1 < p < 00. The Dirichlet problem (Dir-Ak,p) fails to be well 
posed in the sense if k > tan(^), where 1/q = 1 — 1/p. The regularity problem 

(Reg-Ak,p) fails to be well posed in the sense if k < — tan(^). The Neumann 

problem (Neu-Ak,p) fails to be well posed in the sense if k > tan(^). 

In this paper, we demonstrate that one must be careful in specifying in what sense 
well posedness is meant, when considering EVP's for non symmetric coefficients. 
(The notion of well posedness in the sense is defined below.) Indeed, Theorem ll.2l 
below shows that these EVP's can be well posed in an Loo{Lp) sense, as defined 
below, without being well posed in the mentioned sense. 
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To explain these results, we first need to introduce the notion of solutions and 
Loo(yLp) solutions to EVP's. We consider a given divergence form equation 

(1) divA(x)Vf/(t,x) = 

in the upper half plane := {{t,x) G ; t > 0}, where A = (aij)ij=o,i G 
Loo(R; >C(C^)) is a t-independent, complex and accretive coefficient matrix such 
that Re{A{x)v,v) > K|fp, x G R, f G C^, for some k > 0, and where U satisfies 
one of the following prescribed boundary conditions. 

• The Dirichlet problem (DiT-A,p): U{0,-) = u, for a given function u G 
Lp(R; C). 

• The (Dirichlet) regularity problem (Reg-A,p): diU{0,-) = u', for a given 
function u G Wp{K; C). 

• The Neumann problem (Neu-A,^): aoot^of^(0, ■) + (^oidiU^O, ■) = 0, for a 
given function G Lp(R; C). This means that the conormal derivative of U 
is prescribed. 

Throughout this paper, p denotes a fixed exponent such that 1 < p < oo, and q 
is the dual exponent. The regularity and Neumann problems can be thought of as 
EVP's for the gradient vector field F{t,x) = FqCq + FiCi := VU(t,x), rather than 
U itself. Here Cq denotes the vertical basis vector along the t = Xo-axis, and ei is 
the horizontal basis vector along the x = Xi-axis. 

There are two classical methods for constructing a solution U: Lax-Milgram's 
lemma and boundary equation methods. 

SOLUTIONS 

To solve a Dirichlet problem for a sufficiently smooth and localised boundary 
function u{x), one first constructs a function Ui in R^ such that Ui{0,x) = u{x). 
Next, one uses Lax-Milgram's lemma to find a function U2 G ifo(R^), which decays 
at infinity, such that B{U2,ip) = C-ijp) for all ip G i/o(R+), where 

B{U2,i)) ■= jj ^ {A{x)VU2{t,x),Vi){t,x)) dtdx 

and the given functional is £(■?/') := //j^2 {A{x)VUi{t^x)^Vip{t^x)) dtdx. Details of 

this construction for the unbounded domain R^ are found in [H Lemma 1.1]. The 
function U := Ui — U2 now solves equation (P) and has boundary trace u, in a weak 
sense. 

We say that the EVP {Dii-A,p) is well posed in the sense, if for all sufficiently 
smooth and localised u, the solution U constructed above has quantitative bounds 

(2) \mU)\\LAn)<C4u\\L,iny 

Similarly, {Reg-A,p) is well posed in the sense, if for all sufficiently smooth and 
localised u, the solution U constructed above has quantitative bounds 

||iv*(vt/)|U,(R)<c,||«'|U,(R). 

HereiV,(f/)(xo) := sup|^_^,^|<j x)| and iV,(F)(a;o) := sup|^_^^|<i t-iF||L,(Q(t,^.)), 
where Q{t, x) denotes the square centered at {t, x) with sidelength t, are the standard 
(modified) non-tangential maximal functions, and Cp denotes a constant indepen- 
dent of u. 
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Turning to the Neumann problem, this is solved for a sufficiently smooth and 
localised boundary function with J = 0, by applying the Lax-Milgram lemma 
to obtain U e H\Rl), such that B{U,iIj) = for all ip e i^^(R+), where 

i{ip) = — (f){x)ip{x) dx. 
Jk 

Details of this construction for the unbounded domain are found in [H Lemma 
1.2]. This function U solves equation ([T]) and has conormal derivative cj) at the 
boundary, in a weak sense. 

We say that the BVP (Neu-A,p) is well posed in the sense, if for all suffi- 
ciently smooth and localised with J = 0, the solution U constructed above has 
quantitative bounds 

l|iv,(vt/)|U,(R)<Cp||0|U,(R). 

LociLp) SOLUTIONS 

A different method for constructing a solution U to one of the EVP's above is 
the boundary equation method. We are given a kernel function K{t,x;y), which 
for each y G R satisfies the equation ^ in the variable (t,x) G R^. From this we 
obtain, for each auxiliary function h{y) on the boundary, a function 

U{t,x):= / K{t,x;y)h{y)dy 
Jn. 

solving the equation ([T]) in R^. Taking the appropriate trace of U, depending on 
which boundary condition U is supposed to satisfy, we get an equation g = T{h), 
where g denotes either u, u' or 0. If the operator T : Lp(R) Lp{Ti) is an 
isomorphism, then we can solve the equation for h and from this construct a solution 
U. 

Obviously there is a freedom of choice for the kernel function K(t,x;y). In this 
paper, we shall use the boundary equation method from Auscher, Axelsson and 
Hofmann [2]. The Cauchy integral method used here for the Neumann and regularity 
problems actually uses a vector valued kernel K(t, x; y), and constructs the gradient 
vector field F = VU rather than U itself. With some abuse of notation (as the 
invertibility of the boundary equation may depend on the choice of kernel K) , we 
shall say that the EVP's are well posed in the Loo{Lp) sense, referring to the norm 
suPf>o \\U{t, ■)\\p for solutions, if this method gives rise to an Lp invertible boundary 
equation. 

In this paper we shall prove the following surprising, in view of Theorem 11.11 
result. 

Theorem 1.2. The boundary equation method of [2] yields the following result for 
coefficients A^. 

The Dirichlet problem (Dir-Ak,p) is well posed in the Loo{Lp) sense ifk ^ tan(^), 
where 1/q = 1 — 1/p. In this case, the solution Ut{x) = U (t, x) has bound || A^*(?7) \\p < 
C||m||p and convergence \\Ut — u\\p vjhen t ^ 0+. 

The regularity problem (Reg-Ak,p) is well posed in the L^{Lp) sense if k ^ 
— tan(^). In this case, the solution Ft{x) = F{t,x) has bound ||A^*(F)||p < C||m'||p 
and convergence \\Ft — f\\p — > when t — >■ 0^, where fi = u' . 



4 



ANDREAS AXELSSON 



The Neumann problem (Neu-Af:,p) is well posed in the L^oiLp) sense if k ^ 
tan(^). In this case, the solution Ft{x) = F(t,x) has bound ||A^*(F)||j, < C||0||p 
and convergence \\Ft — f\\p when t 0^, where /o + ksgn{x)fi = 0. 

The reason for these seemingly contradictory results is that the two methods con- 
struct different solutions, for some k. To illustrate this phenomenon, we study in 
detail the solutions to the Dirichlet problem in section [H We derive the follow- 
ing explicit expression for the solution to the Dirichlet problem with the boundary 
equation method. 

1 r 2xty + |y|-"Im {(|x| + ttr+^ {y^ - (|x| - itf) } 
3 U{t,x) = - ' ; -u{y)dy, 

where {t^x) G R^. Here tan(7rQ;/2) = k, and a G (1/g — 2, 1/g) is the branch 
obtained with the boundary equation method. Denote the Poisson kernel in ([3]) by 
Pait, x; y). Below the harmonic measures Pq,(0.5, 1; ■) are plotted for some values of 
a. On the other hand, for sufficiently smooth and localised boundary data, the 



k=0 (a=0) k=0.5 (a=0.295) k=2 (a=0.705) 




-4 -2 2 4 6 -4 -2 2 4 6 -4 -2 2 4 6 




solution to the Dirichlet problem is given by (I3|), but with the branch a G (—1,1). 
For k > tan(^), the solution thus differs from the Loo{Lp) solution, and we note 
the following. 
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• The solution uses the branch a e (l/q', 1). Here the kernel Pa(t,x;y) is 
always positive, but does not satisfy the reverse Holder estimates Bg equiv- 
alent with the estimate ([2]), since Pa{t,x]y) ~ ||/|~" around y = 0. The 
solution does not belong to Lp(Iij;) for any fixed t > 0, because of the slow 
decay Pait,x;y) ~ when x ^ oo. 

• The Loo{Lp) solution uses the branch a G (—2 + 1/g, —1). Here the kernel 
Pait, x; y) is not always positive, but does satisfy the reverse Holder estimates 
Bq since also Pait, x; y) ~ when y oo. The solution is not up to 
the boundary in a neighborhood of the origin, not even for smooth boundary 
data. 

The reason why a signed harmonic measure is possible for a G (—2 + — 1), 
without contradicting the maximum principle, is that the solution, even for < 
u e C^(R), satisfies limt^Q+ U{t,0) = —oo. This prevents the L^o approximation 
needed for applying the maximum principle. Indeed, ^ shows that 

N cos(7ra/2) f t^+"\y\-'' , , , 

We summarise the main point of this paper. When considering EVP's with non 
symmetric coefficients, it is important to specify which solution is meant. When 
the solution to the Dirichlet problem does not satisfy ((21), there can still exist 
another solution f/t(x) = U(t,x), which have bounds ||A^*(f/)||p < oo and Lp trace 
\\Ut-u\\p^O,t^O. 

Acknowledgments. The author thanks Pascal Auscher and Steve Hofmann for 
many interesting discussions on the topic of this paper. 



2. Computation of Cauchy integrals 



In this section we explicitly calculate the basic operators we need in order to solve 
the EVP's with the boundary equation method from [2]. As in [21 equation (1.5)1, 
we rewrite the equation ([T]) for U as the equivalent first order system divA^F = 
and curlF = for the vector field F = Vf/. Solving for the vertical derivative, this 
first order system reads 



(5) 



dtF + 



k{sgn{x)d^ - dr,sgn{x)) 

-d.. 



0. 



Throughout this paper we shall identify /oCo + Ziei = [/o, /i]*. The tangential matrix 
operator in ((51) will be denoted T^, and is seen to be a self-adjoint operator in L2(R), 
with domain 

D(Tfc) := {/oeo + /iCi ; /o G H\K), h - ksgn{x)U G H\R)}. 

To solve the EVP's, we need to calculate certain operators in the functional cal- 
culus of the self-adjoint operator T^, in particular we need the following result. 

Theorem 2.1. The Cauchy (singular) integral operators for are 



sgn(Tfc)/(x) 



1 

TT 



—p. V.J ^^z^dy 
p.v. f ^^^dy 



TT 1 + A;2 



r My)+ksgn(y)fi{y) i 

J \x\+\y\ 



sgn(^a; 



kfo(y)-sgn{y)fi{y) 

H+\y\ 



dy 
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and 



e-'i^'=ix+(r.)/(x) = - 



tfo{y)-i^-y)fi(y) 

fi + (x 



th(y)+{^~y)fo(y) 



dy 
dy 



+ 



k 



27r 1 + A;2 



3gn(a;) / 



r t(kfoiy)-sgn(y)My))+{\x\ + \y\){fo{y)+ksgn{y)hiy)) 
J t2+(|^| + |j,|)2 

tifoiy)+ksgn{y)fi{y))+{\x\ + \y\){-kfo{y)+sgn{y)fi(y)) 
i2+(|x| + |y|)2 



dy 



dy 



Here x±{^) denotes the characteristic function of the right/left complex half plane. 
We write sgn(z) := X+i^) ^X-i^) ^ind \z\ := zsgn{z), for z G C. Note that \z\ does 
not denote absolute value for non real z. 

Lemma 2.2. For non real i\ G C, the resolvent {iX — Tk)~^f = u is given by 



sgn(A(x - y))fo{y) - ifi{y))dy 

, J e-M{ - tfoiy) - sgn{\y)My))dy 
_sgn{Xx) J e-^^y^{- ifoiy) - sgn{Xy) fi{y)) dy 



Uo{x) 


SgnA 


- J ^-\X(x-y)\ / 




2 


J ^-\X{x-y)\) 



+ e 



\Xx\ 



2(l-ifcsgn(A)) 



To prove the lemma, we need to solve {i\ — Tk)u = f for m. Thus we are looking 
for u such that 



U'n 



u[ = iXuo- /o 



-iXui + fi, 
fo, 

for X ^ 0, and where u G D(Tfc), i.e. Uq is continuous at x = 0, whereas 

Mi(0+) - ui(O-) = 2kuo{0). 
Multiplying the system of equations with M = [—i, 1; 1, —i] gives the diagonal system 

-Xvo + go, 
v[ = Xvi + gi, 

for V = Mu and g := [—1, — z; i, 1]/. Integrating these equations and using the jump 
condition at a; = gives the formula in the lemma. 

Lemma 2.3. The operators Pt = {l + f^T^)'^ and Qt = tTk{l + t^T^)-\ for t > 0, 
are 



Ptfix) 
Qtf{x) 



'"-^l/*/o(2/) + 



l+fc2' 



<\A+\y\)/t 



( - ^/o(y) + sgn{y)fi{y)))dy 



J J.(^eH-^l/Vi(y) + ^e-(l-'l+l^l)/*(/o(2/) + ksgn{y)My))yy 

I -^sgn(a; - y)e-\''~y\'* fi{y)dy 
J ^sgn(x - y)e-^''-y^/^fo{y)dy 



J 2t H 



-(\A+\y\)/t 



[fo{y) + ksgn{y)fi{y))dy 



I ^^e-(l^l+l^l)/*(fc/o(i/) - sgn{y)f,{y))dy_ 

This follows from Lemma 12.21 and the formulae 

Pt = MiTl-Tk)-'-i^,-n)-') and 
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We are now in position to prove Theorem I2.1[ As in [21 Section 2.3], we use the 
Dunford functional calculus formula 

KTk) = ^. [ &(A)(A-T,)-MA, 
27r« 

where 7 is the boundary of a double sector around R\ {0}. Using b{z) = sgn(z) and 
b{z) = e~*l^lx+(-2) respectively, and Lemma [2?2l gives the formulae in Theorem I2.1[ 
However, the computations can be somewhat simplified by choosing a degenerate 
contour of integration along the imaginary axis. In this case the Dunford formulae 
become 

2 ds 
sgn(Tfc) = - / Qs — , 
^ Jo s 

e-*l^'=lx+m) = - / (Qs cosit/s) + sin(t/.) ) -. 
Changing the order of integration for s and y here, and using that 

/ I p-x/s it/sds _ X ~\~ 
Jo x2 + t2' 

gives the desired formulae. Some of the above computations are of course formal. 
However, they can be justified for example with arguments as in [5]. 

3. Solvability of boundary equations 

In this section we use the Cauchy integrals from Theorem 12.11 to solve EVP's, 
following the boundary equation method described in p]. 

Definition 3.1. Let E'^h := x±{Tk)h be the Hardy projections, with the associ- 
ated Cauchy singular integral operator E^h := sgn(Tk)h. Let the Cauchy extension 
operators be {C^h){t,x) := {e^^\^^\E^h){x), ±t > 0. 

Note that E^ = E^ — Ef^ and conversely E^ = ^(1 ± E^). Given a function 
/i : R on the boundary, applying the Cauchy extension gives a vector 

field F{t,x) = C^h{x) in R^. This is our ansatz for the regularity and Neumann 
problems. Indeed, the vertical derivative of F = e~*l^*l£'^/i is 

dtF = -|Tfc|e-*l'^'=lx+(Tfc)/i = -Tfc(e-*l^'=lx+(Tfe)/i) = -T^K 

Thus F satisfies (JS]), or equivalently the first order system divA^F = and curlF = 

0. This means that F is a gradient vector field F = Vf/, with potential U that 
solves ([T]). 

On the other hand, to solve the Dirichlet problem, we make use of the fact that, 
due to the t-independence of the coefficients Ak, we have 

= dt{divAk{x)F) = divAk{x){dtF) = divAk{x)VFo, 

1. e. the normal component Fq = Cq ■ C^h of the Cauchy extension satisfies the 
equation ([T]). This will be our ansatz for the Dirichlet problem. 

Lemma 3.2. Let h G Lp(R;C2). Then we have bounds \\Ekh\\p < C\\h\\p and 
\\N^:{C^h)\\p < C\\h\\p, and convergence 

\\C+h{t,-)-E+h\\p^O, t^0+. 
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Proof. By breaking up the integrals = + for the second terms in Theo- 
rem ETT], we see that we can write each component of C^h as 



X+ix)[^^ * hi{x) + * h2ix)] +x-ix) 



* h4{x) 



for some Lp functions h^. Similarly E^h is expressed in terms of the Hilbert trans- 
form. The lemma now follows from well known Lp bounds and convergence for these 
convolution operators. □ 

Definition 3.3. Let the double layer potential type operator for Af^ be 



Kf{x) := sgn(x)— p.v. / 

J 

acting boundedly in -Lp(R; C). 



X — y n 



fjy) 

\x\ + \y 



■dy, 



Proposition 3.4. The boundary equation method from [2] constructs solutions to 
the EVP's as follows. 

• A solution to the Neumann problem (Neu-Ak,p) is given by 

VUit,x)=C^[i; 0]\ 

where (p = |(/ + kK)ip. 

• A solution to the regularity problem (Reg-Ak,p) is given by 

VU{t, x) = C+ [-ksgn{x)'ilj{x) ^{x)] * , 

where sgn{x)u' (x) = ^{I — kK){sgn{y)ip{y)). 

• A solution to the Dirichlet problem (Dir-Ak,p) is given by 

U{t,x) = Co ■ [ipix] ksgn{x)ilj{x)y , 
where u = |(/ + kK)*4). 

Proof. Lemma [32] shows that the Cauchy extension has trace limj^o+ C^h = ^{h + 
Ekh). To solve the Neumann problem, we write h = [ijj O]* for the given ansatz. 
Using the first formula in Theorem 12.11 for Ek, the Neumann boundary condition 
becomes 



Co - {Akl{h + Ekh)) 
■■ i [l A;sgn(x 



1 

+ - 
vr 



x—y 



1 k 

nl + k^ 



J \x\+\yry 

sgn(x) / -!^dy 



\^\+\y\ 
U^ + kK^), 



as stated. Similar calculations for the Dirichlet boundary conditions give the stated 
boundary equations. □ 

Remark 3.5. For motivation of the choices made for the ansatzes, we refer to [2]. 
In the language used there, we here express N^E^N^, N^E^N^ and N^E^N^ re- 
spectively in terms of the operator K. The idea is to compress the Cauchy integral to 
a suitable subspace in Lp(R), which is complementary to the null space of the bound- 
ary condition under consideration. We use the projections = [1, /csgn(x); 0, 0], 
= [0, — fcsgn(x); 0, 1], = [1, 0; fcsgn(x), 0], which project onto the complemen- 
tary subspace, along the null space. 
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To prove Theorem II .21 we see that it suffices to investigate the Lp spectrum of K. 
We first note that K = Kq (B Kq m the sphtting -Lp(R) = Lp(R_,_) © Lp(R_), where 
Kq is the operator in the following proposition. 



Proposition 3.6. The operator 



Kaf{,x) := -p.v. / — ::f{y)dy 



77" Jq — iP'' 

is bounded on Lp(R+; C) if and only if a G {—l/p,2 — l/p). Furthermore, if a ^ 
(— 2 — l/p) and k = tan(^a), then I — kK^ is invertible in Lp with inverse 

{I-kKo)-' = j^iI + kK^). 

Proof. We use the isometry 

U : Lp(R+) LpiR) : /(x) ^ e'/^fie'). 

A calculation shows that UKaU^^ = K^+i/p for the convolution operator 

2 e^^ 
TT e^* — 1 



K^g := -p.v.^5^^ * g. 



By standard singular integral theory K-y is bounded, on all Lp spaces, if and only if 
7 e (0,2). This proves the boundedness result for Ka. 
To verify the inverse relation, we need to show that 

{I-kK,/p){I + kK^+yp) = l + k^. 

Applying the Fourier transform, this amounts to 

1 - ki-—- 1 + ki-- — — = 1 + A;^ 
1- zj \ 1 - ;ze*™/ 

where z := e^^^^'^^'^\ This is verified using the relation e*™ = (1 + ik)/{\ — ik). □ 

Proof of Theorem \l.^ Proposition 13.61 shows in particular that I — kK^ is invertible 
in Ivp(R) if 7^ — tan(^). Therefore the boundary equations derived in Propo- 
sition 13.41 are invertible under the hypotheses in Theorem II. 2[ The bounds and 
convergence of the solutions follow from Lemma [X2l □ 

4. The harmonic measure 

In Proposition 13.41 the kernels of the solution operators can now be calculated. 
We here only calculate the harmonic measure Pait,x; ■)dy, where Pa denotes the 
Poisson kernel, i.e. the kernel of the solution operator for the Dirichlet problem. 
The regularity and Neumann problems can be further studied in much the same 
way. According to Proposition 13.41 and Theorem 12.11 a solution to the Dirichlet 
problem is 

TTf. X 1 f f t - {x - y)ksgn{y) \x\ + \y\ \ , ^ 

with 

2 / 2 f°° y^^'^z"" \ 

(6) v(±y) = y^i^y) - ^-p-v- y2 _ ^2 y 1/ > o, 

where k = tan(|a) and a G — 2, 1/g), according to Proposition 13.61 
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Composing these two operators yields the formula ([3]) for the harmonic measure. 
The calculation make use of the following residue calculus formula. 

Lemma 4.1. For a G (—2, 1) and t, ^ > 0, write 
Then, with arg(x + it) G (0, vr), we have 

Z L'- -\- [X — Z j" 

[l-ikY{P -i-l){x + itf+'^ 



-1-ft- /3{x- z) k'^ + l-ft- (3{x + z) 
/«,^,^^r, X, z) - 2 t^ + (x- zY ~ t^ + (x + zY 



Re 



[x + itY — z 



2 



From section [3] it is clear that u \^ ip \^ U gives a solution to (Dir-y4fc,p) when 
a G (1/g — 2, 1/g), and that Ut ^ u rn. Lp(R) when t — >• 0"*". We now investigate 
the Poisson integral formula (JS]) for more general a. First note that Pait, x; ■) is not 
locally Li if a > 1, and if a < — 2 then it does not decay at oo. For these reasons, 
we only consider a G (—2, 1). We shall now consider the branch a G (—1, 1) in ([3]) 
and show that this gives the solution to (Dir-y4fc,p). For the rest of the section, 
we assume that u G C^(R) and a G (—1, 1). 

From (IHl), it is seen that ij G H^/^{K) if M G C^{K). Moreover, we shall use 
the fact that ip := ip — ip{0)e~^^^ satisfies tjj G if^/^(R) and has bounds \iIj{x)\ < 
C min(|x|, l/|a;|)''' for some 7 > 0. 

Lemma 4.2. If u e C^(R) and a G (-1, 1), then 

\VU{t,x)\'^dtdx < 00. 

Proof. Write h := [^{x) ksgn{x)ip{x)y for the ansatz in Proposition 13.41 and note 
that 

Vf/ = V(eo ■ C+h) = dtC+h = -TkC^h = -E+|Tfc|e-*l^'=l/i. 
Since is self-adjoint, it satisfies L2 quadratic estimates. Therefore 

// \VU{t,x)\^dtdx= I ||E+|Tfc|e-*l^'=l/i||^rft 
J Jnl Jo 

rIf 

- [ II m\ni + t'Ti)-\\T,\''^h)\\fi = r^^^~ 

since |2;|^/^e~l^l and + z'^)~^ decays at and 00. Lemma [2731 shows that 

\\h - Pth\\l = \\tp -pt* i^Wl + k'^\\sgn{x)ij -pt* (sgn(x)V^) - 2t{tp,pt)sgn{x)pt\\l, 

where Pt{x) := ^e"'"^'/*. We now recall that ip = ip + 2-?/'(0)pi, where ip G H^^'^{lVj 
and |'?/'(a;)| < C min(|x|, for some 7 > 0, and that 

JR JR F ^ y\ Jn 
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Using the first expression for the norm, we verify that sgn{x)tp E if^/^(R), and 
applying Plancherel's theorem and using the second expression, we show that 

ll/-P**/ll2f <C||/|liv.(R)- 

Thus it remains to show that 

dt f°° dt 

l(V',pOPy + lV'(0)Py \\sgn{x)pi-pt*{sgii{x)pi)-2t{pi,pt)sgn{x)pt\\'^— < oo. 

Here the first term is finite since \ {ip,Pt)\ < C min(t, 1/ty, and an explicit calculation 
for the second term gives || • • • ||^ < C min(t^, 1), which shows that the second term 
is finite. This proves the lemma. □ 

Restricting (l3|) to x = gives the formula ((41). In the quadrants t > 0, ±x > 0, 
the equation divA^Vt/ = reduces to the Laplace equation. For example, for the 
first quadrant, the Poisson integral for the Laplace equation yields 

^^'^^^ = Uo AxH^ + {T-t--y^y ''^y^'^y 

1 f°° 4xts 

when t > 0, X > 0. Inserting the expression 1^ in the second term and calculating 
the integral 



-ds 



7T Jo {4xH^ + -t^ + s2)2) (^2 + ^2) 



+ — Re (l-ik) 



/2 



2 cos(7ra/2) 4x2^2 + (a;2 _ ^2 _ ^2)2 4^^-^ ^ ^ + y 

is seen to give back (l3|). 

It is clear from ((71) that U satisfy the Laplace equation in the quadrants t > 0, 
±x > 0. Furthermore, calculating dxU(t,0±) from ([3]) and dtU(t,0) from (Jll), shows 
that dxU{t, 0+)-dxU{t, 0-) = 2kdtU{t, 0) for t > 0. This proves that divA^Vf/ = 
in R^. Furthermore it is clear from ([7]) that U has boundary trace u in the weak 
sense. This proves that (JS]) gives the solution to the Dirichlet problem when 
a G (—1, 1). Finally we note that it follows from (l7|) and ^ that Pa(t, x;y) >0 for 
all t > 0, x,y E H, when a G (—1, 1). 
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